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> : ABSTRACT 

, We present multiple epochs of Ha spectroscopy for 47 members of the open cluster NGC 3766 

T-H ■ to investigate the long term variability of its Be stars. Sixteen of the stars in this sample are 

I Be stars, including one new discovery. Of these, we observe an unprecedented 11 Be stars that 

■ undergo disk appearances and/or near disappearances in our Ha spectra, making this the most 
^-H ' variable population of Be stars known to date. NGC 3766 is therefore an excellent location to 

study the formation mechanism of Be star disks. From blue optical spectra of 38 cluster members 
and existing Stromgren photometry of the cluster, we also measure rotational velocities, effective 
temperatures, and polar surface gravities to investigate the physical and evolutionary factors 
that may contribute to the Be phenomenon. Our analysis also provides improvements to the 

■ reddening and distance of NGC 3766, and we find EiB ~ V) ^ 0.22 ± 0.03 and {V - My)o = 
I 11.6 ± 0.2, respectively. The Be stars are not associated with a particular stage of main-sequence 

evolution, but they are a population of rapidly rotating stars with a velocity distribution generally 
consistent with rotation at 70 — 80% of the critical velocity, although systematic effects probably 
underestimate the true rotational velocities so that the rotation is much closer to critical. Our 
measurements of the changing disk sizes are consistent with the idea that transitory, nonradial 
pulsations contribute to the formation of these highly variable disks. 

Subject headings: stars: emission-line. Be — open clusters and associations: individual (NGC 3766) 
1. Introduction 
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^Visiting Astronomer, Cerro Tololo Inter- American Ob- 
servatory. CTIO is operated by AURA, Inc. under contract NGC 3766 is a rich, young open cluster in 

the Carina spiral arm that is well known for its 
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many previous studies of this cluster have fo- 
cused on the characteristics of these stars to iden- 
tify their evolutionary status. The cluster has 
bee n the target of numerous photometric stud- 
ies (lAhme dlll96 2t lYilmazlll976 !: 'Shobbrook"l985. 



1987: iMoitinho et al. 1997: iPiatti ct al. 1998; 



Tadrosd[200lt iMcSwain fc Gied l2005bh . But de 



spite these intensive investigations, the cluster's 
age and distance remain somewhat uncertain; 
measuremen ts of its age ra n ge from 14.5 to 25 Myr 
(WEBDAFI: iLvngal Il987t iMoitinho et"all Il997 : 
Tadrosd 2001 ). and its distance is between 1.5 



and 2.2 kpc. The reddening E( B-V) is between 
0.16 and 0.22 (see the discussion of lMoitinho et al "* 



19971) . 



Spectroscopic investigations of NGC 3766 have 
targeted a limited sample of cluster members, fo- 
cusing prima rily on the Be star and superg iant 



cusmg prima rily on tne pe star and superg iant 
populations ( Harris _ 1976 ;_ Merrnilliodl 19821 and 
refere nces therein: ISlettebakI 19851 : iLevesaue et"ar 



20051 ). Even the eclipsing double-lined spec- 



troscopic binary BE Centauri (= HD 100915), 
a member of NGC 3766, has been largely ne- 
glected by rnodern spectroscopic observations 
([Clausen et al.l 120071 and references therein). Eor 
most cluster members, no detailed information 
about their physical characteristics such as tem- 
perature, gravity, rotation, and metallicity are 
known. 

In this work, we present red and blue optical 
spectra for both normal B-type and Be stars in 
the cluster. Like many prior studies of NGC 3766, 
our primary goal is to investigate the Be star pop- 
ulation; but unlike other works, we achieve a more 
complete understanding of this subset of B stars 
by comparing these emission-line objects to their 
non-emission counterparts. Therefore we present 
measurements of the effective temperature, Tcs, 
surface gravity, logg, and in most cases the pro- 
jected rotational velocity, Vs'mi, for 26 normal B 
stars and 16 Be stars in NGC 3766. We use these 
results to improve the known reddening and dis- 
tance to the cluster. Erom multiple epochs of Ha 
spectroscopy, we also investigate the variability of 
the circumstellar disks and estimate the disk mass 
loss/gain rates for 11 Be stars. Einally, we use 



*The WEBDA database is maintained by 
E. Paunzen and is available online at 
|http: / / www.univie. ac. at/ webda / navigation, html | 



the observed disk masses and angular momenta 
to show that nonradial pulsations are a possible 
origin for the disks, and they probably fill dur- 
ing short-lived bursts of mass flow from the stellar 
surface. 

2. Observations 

We obtained spectra of NGC 3766 during mul- 
tiple observing runs in 2003 March, 2005 Eebru- 
ary, 2006 May, and 2007 January— July using the 
CTIO Blanco 4-m telescope with the Hydra multi- 
fiber spectrograph and the CTIO 1.5-m telescope 
with the Cassegrain spectrograph, operated by the 
SMARTS Consortium. The details of all runs are 
summarized in Table [TJ Most of the runs targeted 
the Ha emission line profile with low spectral reso- 
lution to characterize the Be stars' emission; how- 
ever, during one run we observed the blue opti- 
cal region with higher resolving power to observe 
numerous other H Balmer and He I line profiles 
and measure the physical parameters of the clus- 
ter members. 

We selected the targets for each run by giving 
highest priority to the known Be stars in this clus- 
ter (save WEBDA No. 232, which was saturated 
in our photometric study). We then selected other 
B-type stars in the cluster by ranking them accord- 
ing to their y-Ha color to preferentially select any 
weak emissio n stars that were not det ected in our 
photometry (jMcSwain fc Gied l2005allbh . All ob- 
servations were performed by M. V. McSwain ex- 
cept the 2007 CTIO 1.5-m runs, which were taken 
in service mode by a SMARTS observer. Eor the 
Hydra observations, we generally began by taking 
short exposures and then parking the fibers used 
for the brightest stars to avoid saturation in the 
longer exposures. Not all of the known Be stars 
could be observed in one fiber configuration, so 
we took three to four exposures each of two con- 
figurations to observe all of the targets. Therefore 
up to eight exposures of each star were obtained 
with the Hydra runs. We also observed a HeN- 
eAr comparison lamp source just before and af- 
ter the set of cluster observations for wavelength 
calibrations. Eor the CTIO 1.5-m observations, 
we alternated each stellar observation with a Ne 
comparison lamp spectrum. 

The CTIO 1.5-m spectra were reduced and rec- 
tified to a unit continuum using standard routines 
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for slit spectra in IRAF. All of the Hydra spec- 
tra were zero corrected using standard routines in 
IRAF, and they were flat fielded, wavelength cali- 
brated, and sky subtracted in IRAF using the do- 
hydra routine. In comparing the slit spectra and 
fiber spectra for many of the same objects in our 
data set, we find no evidence of systematic differ- 
ences in the background subtraction due to cluster 
nebulosity. For each Hydra spectral configuration, 
we transformed the observations to a common he- 
liocentric wavelength grid and co-added them to 
achieve good S/N for each star. 

The complete sample of stars presented in this 
work is listed in Table [21 Column 1 gives each 
star's ident ification number based on the assigned 
number in iMcSwain k Gie^ (|2005d) : the corre- 
sponding WEBDA numbers are given in column 

2, where available. We obtained Ha spectra for 
each of these stars during at least three epochs in 
most cases, and these are shown in Figures [THSl 

3. Physical Parameters from Spectral 
Models 

3.1. Vsmi Measurements 

We began our investigation of each star's 
physical parameters by generating a grid of 
synthetic, plane-parallel, local thermodynamic 
equilibrium (LTE) atmosp heric models using 
the Kurucz ATLAS9 code (|Kuruczl [l99l . We 
adopted solar abundances and a microturbu- 
lent velocity of 2 km s~^ for these stars, which 
corresponds to the mean microturbulence ob- 
served among_Jate;^tif£e;^jnBii>^^ 



stars ( Lvubimkov. Rostopchin. fc Lambert 2004[ ). 
Each atmospheric model was then used to cal- 
culate a grid of mode l spectra using SYNSPEC 
(|Lanz fc Hubenvll2003l ). 

For the 38 stars with available blue spectra, we 
made a preliminary estimate of their effective tem- 
perature and gravity, Tes and log g respectively, by 
comparing the observed H7, H(5, He I A4143, and 
He I A447I-|-Mg II A4481 hue profiles to our grid 
of Kurucz spectral models. To measure Vsmi, 
we compared the observed He I line profiles to the 
model profiles convolved with a limb-darkened, ro- 
tational broadening function and a Gaussian in- 
strumental broadening function. We determined 
the best fit over a grid of values, spaced 2 km s~^ 
apart, minimizing the mean square of the devia- 



tions, rms^. The formal error, AVsmi, is the off- 
set from the best-fit value that increases the rms'^ 
by 2.7rms^/iV, where N is the number of wave- 
length points within the fit region. Our measured 
Vsmi and AVs'mi are listed in columns 3-4 of 
Table [3 

Even the He I lines may contain some weak 
emission in Be stars, partially filling and narrow- 
ing their line profiles. Furthermore, a number of 
the Be stars show evidence of narrow "shell" line 
components (formed in the outer disk), and the 
presence of a shell component may make the pro- 
file appear too narrow in some cases. Therefore 
we consider our V sin i measurements for the Be 
stars to be lower limits. However, we note that 
the He I lines do not exhibit obvious signs of emis- 
sion among most of the Be stars, and these lines 
are much less susceptible to emission than the H7 
or H(5 lines. For the case of No. 154, a shell star 
with strong emission and contamination present 
in the He I lines, we used the Mg II A4481 line to 
measure Vsini. 

There are few previous measurements of V sin i 
for members of NGC 3766 in the literature, but we 
found th at 10 stars in our sample were also mea- 
sured by ISlettebakI (Il985l) . Our V smi measure- 
ments generally agree well, with the exception of 
No. 154. Slettebak found Vsini = 220 km s'^ for 
that star, nearly double our measured value. We 
emphasize that the exceptionally strong He and 
metal lines of this shell star amplify the difficulty 
of measuring its Vsini. 

3.2. Tcff and log g Measurements of B stars 

For the B stars with Tcff < 15,000 K, we 
used t he "virtual star" method of Huang fc Gied 
()2006l ) to improve our Toff and log g measurements. 
(Their virtual star is a spherically symmetric star 
with constant Tcs and log 5 across the stellar sur- 
face.) They generated detailed H7 line profiles 
using line-blanketed, LTE Kurucz ATLAS9 and 
SYNSPEC codes. Huang & Gies show that the 
H7 line strength and equivalent width can be used 
as starting parameters in a line profile fit to obtain 
unique values of T^s, log 5, and their correspond- 
ing errors. We used their procedure to measure 
these quantities from our observed H7 line pro- 
files. 

Among the hotter B-type stars, non-LTE ef- 



3 



fects alter the equivalent width of the H7 line, 
and thus the LTE Kurucz model line profiles sys- 
tematically underestimate Tcff- Therefore we used 
the new TLUSTY BSTAR2006 grid of metal Hne- 
blanketed, no n-LTE, plane-parallel, hydrostatic 
model spectra ( Lanz fc Hubenv 20071 ) to measure 
Tes and logg for those stars with T^g > 15, 000 K. 
We used their models with solar metallicity and 
helium abundance and a microturbulent velocity 
of 2 km s~^. The grid includes Tes from 15,000- 
30,000 K in increments of 1,000 K and \ogg from 
1.75-4.75 in increments of 0.25 dex. For these hot 
stars we measured T^g and log g by comparing the 
H7 line profile to the rotationally and instrumen- 
tally broadened model spectral line profiles at each 
value in the grid, minimizing rms^ across the line 
region. We then refined our measurements to a 
higher precision using a linear interpolation be- 
tween the available line profiles in the grid. Fi- 
nally, we determined the errors, ATch and A log 5, 
from the values which produce a rms^ no more 
than 2.7rms^/7V greater than the minimum rms^. 
Our measurements of Toff and log 5, with their cor- 
responding errors, are listed in columns 5-8 of Ta- 
ble H 

Many of the B stars in our sample are rapidly 
rotating with Tasini > 200 km s~^. Such a 
rapidly rotating star will experience strong cen- 
trifugal forces that distort the star into an oblate 
spheroid al shape, as recently found for the star 
Regulus (jMcAlister et al. II2OO5I) . The surface grav- 
ity at the equator can therefore be much lower 
than at the poles. For such a rapid rotator, 
our measured Tcff and logg represent the aver- 
age across the visible stellar hemisphere and are 
therefore biased toward lower values, causing the 
star to appear more evolved. Because the polar 
regions are not distorted, the surface gravity at 
the poles is a be tter indicator o f the e volutionary 
state of the star. Huang fc Gied ( 20061 ) performed 
detailed spectroscopic modeling of such distorted 
rotating stars to determine a statistical correction 
factor for logg, averaged over all possible i, for a 
variety of stellar models. We made a bilinear in- 
terpolation between their models to convert our 
measured logg to loggpoiar for a more accurate 
comparison between slow and rapid rotators. 

For each B star, we also measured its mass. 
Mi,, and radius, i?^, by interpolating between the 
evolutionary tracks for non-rotating stars from 



Schaller et al] (|l992l ). The errors AM* and AR^ 
correspond to our measured ATeff and Alogg. 
Our results for loggpoiar, Af*, AM*, i?*, and AiZ* 
are also listed in Table [31 columns 9-13. 

3.3. Toff and log 5 Measurements of Be 
stars 

Star No. 196 did not show any sign of Be emis- 
sion until our most recent observations, hence we 
included it among the normal B stars and mea- 
sured Toff and logg from the H7 line in its 2006 
blue spectrum. For other Be stars in our blue spec- 
tra, the above method to measure Teff and logg 
was not useful because the equivalent width of the 
H7 line may be decreased by emission in the line, 
even if the line profile does not exhibit obvious 
signs of emission. 

Our first attempt to measure the Be stars' Toff 
and logg reUed upon the He I AA 4143, 4388, 
4471 lines. We fit these three hnes using TLUSTY 
model spectra using the same procedure described 
above for the H7 line. However, these He I line 
strengths are less sensitive to Toff, and the line 
wings show only a very small dependence on logg. 
Therefore these line fits resulted in very large er- 
rors for both parameters in many cases. Further- 
more, the Be disks contribute continuum flux that 
dilutes the apparent strength of the He I lines, so 
logg values measured from these lines do not al- 
ways agree well. 

To improve our Toff measurements for the Be 
stars, we turned to available Str6m£ren_mi_;_Ci , 
and indices for our targets (|Shobbrookl Il985l 
1987t WEB DA). Several temperature relations 
for Stromgren indices are available in the liter- 



ature (see 'Napiwotz ki et al.l Il993l and references 
therein) , so we began by using our normal B-type 
stars with Toff > 15000 K as calibrators to deter- 
mine the best relation for our data. We adopted 
the reddening value of E{b — y) = 0.15 (corre- 
sponding to E(B - V) = 0.2) for NGC 3766 
(jShobbrookl Il985l : iMoitinho et all Il997tl . Eight 



stars in our sample (Nos. 16, 42, 49, 54, 57, 161, 
170, and 196) have available Stromgren indices 
as well as Toff measured from our H"f line fits, 
and Napiwotzki et al. provide eight calibrators 
with well-known Tcff > 15000 K (measured from 
their absolute integrated stellar flux) and available 
Stromgren data. For these 16 B stars, we found 
the best overall agreement using the temperature 
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relation from iBalonal (|l984[ ). shown in Figure [H 
However, we found that Teaiona systematically un- 
derestimated Toff, and a correction factor was nec- 
essary to improve their agreement. We performed 
a linear fit to the data and found the relationship 



TcfT,fit = 1.052 TBalona - 359.636 K 



(1) 



(excl uding one outlying point from lNapiwotzki et al 
1993( ). The slope of this correction is vir- 
tually identical t o the valu es found by both 
Gies fc Lambert ( 1992[ ) and Cunha fc Lambert 
(|l994l ). After applying this correction, we found a 
mean scatter of 264 K from our two independent 
measurements of T^ff for the B stars in NGC 3766. 

Tkaiona relics upon the dereddened cq index as 
well as the narrow-band f3 magnitude, so deter- 
mining an accurate temperature for the Be stars 
also requires confidence in p. However, the Be 
stars' emission makes the P magnitude highly un- 
reliable. Therefore we used the B star calibrators 
from our sample (listed above) to inves tigate sev- 
eral ( cp,/?) relations in the liter ature l| Crawford 
19781 : iBalona &: Shobbroo"^ll984h . The (co,/3) di- 
agram is essentially a H-R diagram that reveals 
temperature and evolutionary trends in a popula- 
tion. In Figure El we show that the values of cq 
and /3 for these B-type calibrators generally agree 
with the r elations for luminosity class III and V 
stars from Balona fc ShobbroolJ ( 1984 ). However, 
the calibration stars have 3.43 < \ogg < 3.98 
(3.49 < loggipoiar < 4.17) since these hot stars are 
evolving along the MS, and neither (cq, f3) relation 
can be applied to the entire population. Therefore 
we performed a linear fit to account for the range 
in evolution, and we found the relationship 



= 0.417 Co + 2.545 



(2) 



among the B stars in NGC 3766 with Toff > 
15000 K. The mean scatter between Pat and the 
measured /3 is 0.010, which implies an additional 
error of 144 K in Teff using the corrected Teaiona 
relation above. 

IShobbroo^ ( 1985 . 1987) presented Stromgren 
photometry for nine Be stars that have accom- 
panying blue spectra in this work (excluding No. 
196, which we used as a calibrator) and four Be 
stars with only red spectra in this work. Thus 
we measured Teff for all 13 of these Be stars us- 
ing the adopted Pm and the corrected Teaiona 



as described above. We adopt a total error of 
(2642 -t- 1442)"-5 = 301 K for T^a measured with 
this method, and the results are listed in columns 
5-6 of Tabled 

Using this technique, only the "shell" Be star 
No. 154 results in Toff < 15000 K . The compar- 
isons between Tsaiona and our independently mea- 
sured Tos are less reliable below this temperature 
threshold, hence we used only B star calibrators 
with Tcff > 15000 K. However, spectroscopy is 
even less likely to produce accurate measurements 
for No. 154 since its spectrum is contaminated by 
He I emission and the metal lines are exception- 
ally strong due to the disk's edge-on orientation. 
Therefore we include in Table 5] its Tcff from the 
corrected relations of iBalonal (jl984l ). but the er- 
rors are somewhat higher than for the other Be 
stars. 

The values of log g are more strongly dependent 
on P, so we were reluctant to use our P fg with 
the St romgren relation for log g given by I Balona 
(jl984l ). Instead, we used the calculated T^g to 
determine the Be stars' bolome tric corrections, 
BC, from iLanz fc Hubenvl ( 2007 ). at first assum- 
ing \ogg = 4.0. For the cooler star No. 154, we 
interpola ted t h e BC from the values for MS stars 
given 

by|c3 dimi. We calculated each stellar 

radius, i?*, and luminosity, Lj,j_jising_the mea- 
sured Toff, BC, V magnitude ("S hobbrookl Il985l 
19871 : WEBD A), distance modulus {V - Mv)o 



11.73 ± 0.33 dMoitinho et al.lll997l). and E(B 
V) = 0.2 ± 1 (|Shobbrookl 119851 : iMoitinho et al 



19971 ). We measured the stellar mass, M^, from 



the computed Toff and T^ by int erpolating be- 



tween the evolutionary tracks of iSchaller et al 
(jl992[ ). Finally, we obtained a preliminary value 
of \ogg from and i?^. Since the BC is weakly 
dependent on logj;, we improved the BC from the 
initial estimate and iterated to compute the final 
logg. We adopt a formal error in logg computed 
from AToff and the quoted errors in {V — Mv)o 
and E{B ~ V). 

We note that the scatter between the calibra- 
tors' log^poiar and their Stromgren logg is identi- 
cal to the formal error. Based on this good agree- 
ment, we do not perform any further correction 
to obtain log ^poiar for the Be stars measured with 
this technique. However, we show below that the 
Be stars are more rapidly rotating than the nor- 
mal B stars, which may make them appear slightly 
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more evolved and artificially brightened. Thus 
these values of log gpoiar are lower limits. The final 
parameters for these Be stars are listed in Table 

El 

Two Be stars in our sample of blue spectra have 
no available Stromgren data in the literature. By 
coincidence, these stars (Nos. 31 and 127) have 
two of the best Toff and log g measurements from 
our preliminary He I line fits. Therefore we adopt 
the mean measurements from the He I AA 4143, 
4388, 4471 lines, with a formal error determined 
from the values which produce a rms^ no more 
than 2.7rms^/iV greater than the minimum rms^. 
A contour plot of the errors in Tcff and logg for No. 
31 is shown in Figure [SI and the final values are 
listed in Table [3l We determined their log gpoiar 
using the same method as the B stars. 

3.4. Discussion of Physical Parameters 

The resulting values of Toff and loggpoiar 
are plotted in Figure [7] with the corresponding 
evolutionary t racks for non-rotating stars from 



Schaller et al 



(1992). The symbol sizes are pro- 
portional to Tasini to investigate the relation 
between rotation and the apparent evolutionary 
state of each star, but no trends are observed. 
Also in Figure [3 we p lot the isochrones fo r 25-50 
Myr populations from Leieune &: Schaereil (j2.001). 
Our distribution of T^s and log gpoiar are generally 
consistent with a population within this range, 
indicating an age slightly greater than previous 
estimates for NGC 3766. 

The masses of the rapidly rotating stars may 
be overestimated since the evolutionary tracks do 
not account for rotation. Wc minimize this effect 
by using the corrected log^poiar for all stars, but 
slight mass differences may still be present since 
rapid rotation is expected to alter the evolutio n 
dHcgcr & Langcr '200(1 iMevnet fc Maedeil l2000f) . 
[ilcAhstcr et al. {2003) found a small, 15% mass 
discrepancy for the rapid rotator Regulus when 
comparing its mass derived from non-rotating evo- 
lutionary tracks and its true mass from a detailed 
spectroscopic and interferometric analysis. 

We find that the Be stars are generally among 
the hotter stars in our sample, consistent with our 
earlier findings from a photometric investigation 
of 48 open clusters that Be stars are preferentially 
found among the more luminous cluster members 



()McSwain fc Gied l2005bll . There are no system- 
atic differences between the loggpoiar of the Be 
star and normal B star populations of the clus- 
ter, indicating that the Be stars in NGC 3766 are 
distributed across a range of log (?poiar and are not 
associated with any particular stage of the MS evo- 
lution of B-type sta rs. This is also consistent with 
our earlier result s ( McSwain &: Gied 2005b ) and 
with lZorec et al.l ( 20051 ) . who performed an evolu- 
tionary study of field Be stars . On the other hand, 
Levenhagen fc Leisteil ( 2006 ) found that field Be 
stars are preferentially found at later stages of the 
MS evolution. They did not perform any correc- 
tions for gravity darkening among these rapid ro- 
tators, so their Be stars may be found closer to 
the zero- age MS than their results suggest. 

4. Rotational Velocities 

Be stars are often described as a population of 
rapidly rotating B-type stars with a true rotational 
velocity V^ot comparable to the critical velocity 
Vci-it ■ (Note that V^t should not be confused with 
V in our measured Tasini, which may be subject 
to systematic errors such as gravitational darken- 
ing or weak emission i n the He I line s.) How- 
ever, a recent study bv ICranmei ( 20051 ) has cast 
some doubt on their fast rotation. He compared 
available V sin i measurements in the lYudinl (|200ll ) 
database of Oe, Be, and Ae stars to the predicted 
distribution of V sin i accounting for gravity dark- 
ening, limb darkening, and observational effects. 
He found that the Be stars of the Yudin database 
have intrinsic rotations between 40-100% criti- 
cal, with more early-type Be stars having signif- 
icantly subcritical rotation. However, the Yudin 
database is a compilation of measurements from 
many different authors and instrumental setups, 
and therefore the available V sin i may contain sig- 
nificant systematic differences. Our measurements 
of V sin i rely upon data of identical origin and 
measurements of the same spectral lines in nearly 
every case. Therefore we have greatly reduced the 
systematic differences among our Be star measure- 
ments, and our measurements are also a reliable 
comparison of Be stars relative to normal B stars. 

A rotationally distorted star has an equatorial 
radius i?c — 1.5i?p in the Roche approximation; 
here, Rp is the star's polar radius. For simplicity, 
we assume Rp is equal to the radius i?* derived 
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from the position in Figure [7] and given in Table 
[3l The resulting critical velocity is 



(3) 



and is included in Table [51 A small mass discre 
ancy, as found for Regulus (jMcAlister et al. | [200" 
will not affect V^crit significantly since a 15% error 
in mass produces only a 4% error in Vcrit- Here 
we investigate the rotational properties of several 
edge-on Be stars in our sample and compare the 
distribution of Be star velocities with the normal 
B stars to compare the two populations. 

Two Be stars in our sample, Nos. 92 and 139, 
show Ha emission only in the line wings, with a 
deep absorption profile, suggesting the disks are 
observed nearly edge-on. This is probably a good 
assumption for No. 139 since it has F/T4rit > 0.7. 
However, No. 92 may not be edge on since i ts Ha 



profiles between 1985-1990 (|Balona et al.Ml99lD 



do not resemble the profiles we observe. Either 
the disk is precessing or other structural changes 
have occurred. A third star, No. 154, exhibits a 
shell spectrum with strong metal lines that also 
suggests an edge-on orientation. 

If we assume that all of the Be stars have 
V/Vcrit = 0.95, then our measured Vsmi indi- 
cate that stars 92, 139, and 154 have i = 34°, 
48°, and 15°, respectively. Assuming a slower, 
V/Vcrit = 0.70 increases the derived values to 
i = 49°, 90°, and 20° respectively, but the low i for 
Nos. 92 and 154 are still inconsistent with the ob- 
served line profiles. However, if we fix i = 80° for 
all three stars, we find that y/V^rit = 0.53 for No. 
92, V/Vcrit = 0.71 for No. 139, a nd V /Vcrit = 0.24 



for No. 154. iTownsend et al.l (|2004l ) show that 



measured Vsini may be too low due to gravita- 
tional darkening in highly distorted, rapidly ro- 
tating stars. But even accounting for a 20-33% 
underestimate in Vsini for Nos. 92 and 154, we 
do not find that these two Be stars are near critical 
rotation. 

The distributions of V sin i for both the Be stars 
and normal B stars are ploted in Figure [H Al- 
though this sample of 38 stars is small, this study 
is the first to measure V^sini for both groups us- 
ing a consistent method and a uniform data set. 
Our measured V sin i for the Be stars may be only 
lower limits, but we clearly find that the Be stars 
of NGC 3766 are more rapidly rotating than the 



normal B star population. In Figure [H we also 
include theoretical distributions of two uniform, 
rapidly rotating populations with V — 0.7 Vait 
and V — 0.8 V^^rit, assuming V^^it — 386 km s~^ 
(the mean value among the Be stars, excluding 
No. 154 due to its large measurement errors). 
The distribution of Be stars is consistent with 
0.7 Krit < < 0.8 T4rit in most cases. However, 
if each measured V sin i is underestimated by 20- 
33% due to grav i tation al darkening, as claimed by 
Townsend et al.l (|2004r ). the distribution of Be star 
velocities is consistent with T^ot > 0.84 Vent- Nos. 
92 and 154 represent significant exceptions to this 
rule, as discussed above. 

5. Other Results from Spectra 

We do not wish to broaden the scope of this 
paper to discuss the radial velocities and helium 
abundances observed in our spectra. However, we 
found several instances that suggest these topics 
are worth exploring in future studies of NGC 3766. 

Figure [1] includes the variable Ha profile of 
No. 61. The absorption line is deep and narrow 
in 2005, yet wide and shallow in 2006 and 2007. 
These line profile changes are not consistent with 
a Be star, and they probably indicate line blend- 
ing in a double-lined spectroscopic binary (SB2). 
Likewise, No. 197 may be an SB2. We observed 
double Mg II 4481 hues in 2006, and Figure 3 
shows a large wavelength shift in its Ha line profile 
between 2005—2006. There is possible line blend- 
ing visible in the 2005 and 2007 Ha profiles. 

We measured the physical parameters of the B 
stars using spectral models that have solar helium 
abundances, and most were consistent with a solar 
abundance. However, we noticed a large number 
that have abnormal He abundances; Nos. 41, 55, 
94, 118, and 178 appear to be He strong, while 
Nos. 126, 129, 170, 173, and 197 appear He weak. 
Since all cluster members are expected to have the 
same He abundance, this may indicate widespread 
magnetic fields in the cluster which can produce 
nonuniform distributions of He a cross the stel- 
lar p hotospheres (see discussion of Huang fc GiesI 
20061 ). The He weak stars may also be binaries 
with cooler companions adding fiux but diluting 
the appearance of the He I lines. 

Finally, we used our spectra to reinvestigate 
the reddening and distance to NGC 3766. For 
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the nonemission B stars w ith loggnoi a r > 3.9, 



we used the MS relation of iHarmaned (|l988l ) to 
assign a spectral type to each star based on its 
measured Tes- We then used the MS relation of 
IWegneij (1994) to find the stars' intrinsic colors, 
{B — V)o- Observed B — V colors were obtained 
from WEBDA. For the 17 B-type, MS stars in our 
sample, we found E{B -V) = 0.223 ± 0.030. We 
then used the B stars' V magnitude (WEBDA) , 
BC ()Lanz fc Hubenvll2007l [Malagnini et al.lll986h . 
and our measured i?+ to compute the distance 
modulus, {V - My)o. The resulting {V - My)o = 
11.42 ± 0.15. To determine the mean {V - Mv)o, 
we excluded No. 45 with a somewhat lower {V — 
My)o of 10.3, probably a foreground star. 

While the B stars of NGC 3766 are rotating 
more slowly than the Be stars, they have a mean 
ysini/Vcrit = 0.5. Assuming a random distribu- 
tion of i, this corresponds to a mean y/Vdit = 0.7. 
Such rapidly rotating stars are probably rota- 
tionally distorted with surface areas « 1.2 times 
larger than spherical star with R^, = i?poiar, caus- 
ing them to appear overluminous. The distance 
modulus is then underestimated by a factor of 
2.5 log 1.2 = 0.21. With a realistic distribution 
of V and i, this correction may not be appropri- 
ate for all stars, and the range in A{V — My)o 
is 0.02 - 0.30 for 0.28 < V/Vcrit < 0.83. We ap- 
ply this correction to find a final {V — Mv)o = 
11.6 ± 0.2 for NGC 3766, corresponding to a dis- 
tance of 1.9 — 2.3 kpc. Our new values of E{B~V) 
and {V — My)o are highly consistent with previous 
results for this cluster. 

6. Be Star Variability 

While the long term Ha variability of Be 
stars has often been not ed in the literature (e.g. 
Porter fc Riviniu"i 2003 ). few studies have at- 
tempted long term monitoring of a large sam- 
ple of Be stars to qu antify their variability. 



Hubert fc Floquetl (|l998l ) investigated 273 bright 
Be stars with y < 7.5 that were observed by 
Hipparcos between 1989 August and 1993 Au- 
gust. They found that early type Be stars exhibit 
a very high degree of variability while most late 
Be stars maintained a constant magnitude during 
the duration of the Hipparcos mission. They also 
identified 14 Be stars with recurrent short-lived 
outbursts (0.06 < AHp < 0.3 over timescales 



of 50—500 days) and 8 with long-lived outbursts 
{AHp > 0.12 over timescales of > 500 days) in the 
sample. However, the Hipparcos magnitude. Hp, 
covers a broad waveband ranging from 3400 to 
8500 A, and it traces large scale variations in the 
continuum from scattered light associated with 
the disk. Ha spectroscopy is better suited to mea- 
suring variations over a larger dynamic range in 
the disks of Be stars of all spectral types. 

Figures [T]-[3] show the Ha line profile variations 
of all 47 stars in our spectroscopic sample. We 
were intrigued to discover a new Be star (No. 130) 
in NGC 3766 based on our comparison of its Ha 
line profiles. In 2005 February, the absorption line 
was broad and shallow, with some bumps that sug- 
gest a weak emission disk, nonradial pulsations, or 
simply noise. By 2006 May, however, the depth of 
the line had increased while its width remained 
nearly the same, but a slight asymmetry in the 
line suggested that perhaps a very weak disk may 
still have been present. The line remained in ab- 
sorption during our most recent 2007 observations. 
While an inverse relation between the line width 
and depth is associated with line blending in an 
SB2, the observed variations in No. 130 cannot 
be attributed to such line blending. Instead, the 
three sets of observations indicate the presence of 
a weak circumstellar disk in 2005 that largely dis- 
appeared by 2006. 

In fact, we were startled to find that a total 
of 11 Be stars in our sample exhibit significant 
changes in their disk state between 2003-2007. 
Stars 25, 31, 73, 83, 92, 98, 119, 130, 133, 139, 
and 196 each show Ha in absorption in at least 
one of the four epochs available, yet each star also 
experiences at least one epoch when the Ha line 
is partially or fully filled with emission. These ob- 
servations suggest that the mechanism responsible 
for the disk formation is unstable over timescales 
of only a few years, c onsist ent with the results 



of [Hubert fc Floquetl (|l998f ). The disk growth 



of No. 31 during 2007 reveals that the disk for- 
mation can be very rapid, only requiring a few 
days or weeks, as observ ed in other Be stars (e.g. 



Grundstrom et al.ll2007f l. However, the high frac- 



tion of Be stars that show significant variability 
in a single cluster is unprecedented, and the clus- 
ter should be monitored with increased frequency 
to measure an accurate timescale of the disk state 
changes. 
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Five additional Be stars (47, 127, 154, 198, and 
200) show Ha emission in every observation, al- 
though the emission strength is usually variable. 
Such variations in Ha emission strength are typi- 
cal among Be stars. 

To quantify the observed changes in the Be star 
disks, we measured the equivalent width, Wx, of 
the Ha line in each of our red spectra by nor- 
malizing each spectrum to a unit continuum and 
integrating over the line profile. Columns 3—6 of 
Table [2] give the measured values of W\ for each 
spectrum. The error in each measurement is about 
10% due to noise in the continuum region. 

With multiple epochs of Ha observations avail- 
able for so many variable Be stars, it is worthwhile 
to esti mate the changing size of t he circumstellar 
disks. iGrundstrom fc Gied (|2006h describe a sim- 
ple model to measure the ratio of the projected ef- 
fective disk radius to the stellar radius, i?disk/-R*, 
and the density at the base of the disk, pQ, using 
W\ and Teff as input parameters. We determined 
the disk inclination i by assuming that each star is 
rotating with V at 70% of Vait in most cases. We 
also assumed a disk truncation radius of 100 i?*, 
the nominal value unless a close binary companion 
is present. The disk temperature is assumed to be 
constant at 0.6 T^e- 

To estimate the total masses of the disks, we 
used an axisymmetric, isothermal density distri- 
bution. 



p{r, z) = po[—] exp 



1 



2 \H{r) 



(4) 



( Carciofi fc: Biorkmanll2006h and a radial density 
exponent n = 3 , typical of other Be star disks 
(iGies et al.ll2"007h . The scale height of the disk is 



where 



Hir) = Ho 



Hn 



Vera 



R*, 



I kT 
pmn'' 



(5) 

(6) 
(7) 



The resulting disk properties are listed in Table 
m for each Be star with an available Toff measure- 
ment. The time of each observation is listed in 
column 2, and the dates for the CTIO 4m-|-Hydra 
observations are less precise since these were 
measured from co-added spectra rather than indi- 
vidual spectra obtained with the CTIO 1.5m tele- 
scope. Our determined values for i, poi Rdisk/R*, 
and Mdisk are listed in columns 4-7. 

Certainly, these disk measurements should be 
viewed with caution since the true density profiles 
may be quite different from the assumed distri- 
bution. The disk density exponent may be dif- 
ferent from our assumed value of n = 3, and the 
resulting mass estimates may differ from our re- 
sults by orders of magnitude. Contemporaneous 
optical and infrared spectra have also revealed ev- 
idence of densi ty waves that alter the azimuthal 



and (3 = 1.5 for an isothermal disk ( Biorkman fc Carciofil 
20051 : ICarciofi k BiorknlM]|2006h . 



disk structure (jWisniewski et al.l 120071 ). Finally, 



the disks are likely not isothermal as we assumed, 
and t he true thermal structure m ay be very com- 
plex ( Carciofi fc Biorkman 20061 ). Therefore the 
processes associated with the Ha emission pro- 
file are r nuch more complex than a ssumed in the 
model of Grundstrom fc Giei ( 20061 ). but by using 
their simple model we obtain relative estimates of 
the variations in disk base density that are sensi- 
ble provided the disk density exponent is assumed 
constant in time. 

7. Disk Formation by Nonradial Pulsa- 
tions 

The estimates for po and Mdisk above are not 
sensitive to i and offer a consistent method to 
measure changes in the size of the disks. Thus 
we can estimate the disk growth rate (or dissipa- 
tion rate), AMdisk/Ai. Our generally sparse ob- 
servations cannot determine whether A A/disk/ At 
is constant or highly variable over timescales less 
than one year; however there is an indication from 
the two closely spaced observations of Nos. 25 and 
31 that the disk size can change rapidly. We pro- 
vide our measured AMdisk/Ai between each pair 
of observations in Table ID The long-term buildup 
rates are comparable for every Be star in our sam- 
ple, hinting that all of the disks might be formed 
by the same mechanism. Finally, the buildup and 
dissipation rates are also comparable in magni- 
tude. 
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Many Be stars are known to exhibit nonradial 
puls ations (NRP) with the mode £ = 2, m = 
±2 ( Rivinius et al. 2003f l. and such pulsations are 
commonly proposed as the sour ce of kinetic energy 
to inj ect material into the disk ( Porter fc Riviniua 
l2003h . The m = +2 pulsational mode is retro- 
grade, and naively it will counteract the rotational 
velocity of the star and hinder disk formation. 
However, iTownsendl (|2005) shows that retrograde 
mixed modes behave in a way that could con- 
tribute to Be disk formation. Although their phase 
velocity is retrograde, their group velocity is pro- 
grade. The density enhancements occur when the 
pulsational velocity perturbation is in the same di- 
rection as the rotation - a configuration favorable 
to mass ejection. He also found the mixed mode 
instability strip likely overlaps with the tempera- 
tures and spectral types of known Be stars if their 
rotation is nearly critical. 

Observations of known (3 Cephei pulsators and 
slowly pulsating B stars indicate that their pulsa- 
tion modes have velocity amplitudes on the order 
of 10 km (combining the radial, azimuthal, 
and longit udinal cornpone n ts of the total veloc- 
ity vector; iDe Riddeij[200ll) . lOwockil (|2005h shows 
that the velocity needed for ejection into the disk 
is AV^rb = VcTit — Vrnt- FoT a weak atmospheric 
process such as NRP, the atmospheric sound speed 
a ~ 12 km s~^ must be comparable to AT^ib- If 
we assume that the stars have Vrot = 0.95 Vcru, 
then AV^ib ~ 20 km s^^, highly comparable to 
a and the observed pulsation velocities in other 
NRP stars (jbc Riddcr 2001). As we show above, 
this assumption of near critical rotation is reason- 
able for most of the Be stars in NGC 3766. 

To compare the NRP velocity with the outflow 
velocity required to fill the Be disks, we use the 
equation of mass continuity for an equatorial out- 
flow with velocity t;r,oq: 



cr(r) 



AMdi.k/At 

27rr Vr on 



(8) 



(9) 



Since we remove the vertical dependence of the 
disk density by using its surface density, this 
method effectively assumes that all of the mate- 
rial leaves the star through a cylindrical surface 



where the surface density of the disk is 



of arbitrary height. For a disk reaching down 
to the stellar surface, the mass must flow across 
this surface, so we use the surface density at the 
stellar surface, a{Rcq), in the continuity relation. 
Since the £ = 2, m — ±2 modes have the great- 
est pulsational amplitude at the stellar equator, 
where the effective surface gravity is also a min- 
imum, it is reasonable to expect the majority of 
the disk material to be ejected from the equator 
in the NRP model. The observed disk formation 
rates are generally AMdisk/Ai = lO"^^ - 10~" 
Mq yr~^ with po = 10~^^ g cm~^ from Table [H 
This implies a very slow surface flow on average, 
Vr,cq ^ -01 — 0.1 km s~^. However, a few observa- 
tions reveal that the surface flow can be an order 
of magnitude faster, but this is still consistent with 
NRP. 

Likewise, we use our observed disk densities to 
demonstrate that NRP are a capable source of an- 
gular momentum, L, for disks. The rotational in- 
ertia, /, of a thin disk shell is given by 



I{r) = 27rr'^cr(r) dr 



(10) 



and the total angular momentum of the disk, ex- 
tending from the stellar surface to lOOi?*, is 



L = I{r) u!{r) dr. 



(11) 



^•^ is the angular ve- 
We flnd that L - 
lO''^ g cm^ s~^ for the observed disks, with 



? —9 

g cm s . 



Osakil (|l986h 



where Lu{r) — (GM*)°-^ r 
locity of a Keplerian disk 
10^3 

AL/At - 10^2 - 10^6 
shows that the angular momentum flux in the 
equatorial plane due to NRP with £ = 2, m = ±2 
is given by 

^ = 27ri?2^ (7(i?eq) A^k sin 5. (12) 

Here, A and 2Ak are the radial and azimuthal 
amplitudes of pulsation, respectively, and S is the 
phase shift between the two velocity components. 
Thus we find that a small azimuthal pulsation with 
2Ak ^ 0.01 — 0.1 km s~^ can, in principle, pro- 
vide enough angular momentum to eject material 
into the observed disks, and both the radial and 
azimuthal amplitudes are highly consistent with 
the observed pulsation amplitudes of other B stars 
with NRP (,De Riddeill200ll ). 

Our results suggest that the NRP may be 
a transitory phenomenon for the Be stars, and 
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the disks may fill substantially during short pe- 
riods of surface activity. Evidence for changing 
pulsational modes among s everal stars in NGC 
3766 has been observed by iBalona et al.l (Il99lh. 
van Vuuren et al. I (|l988l ). and lBalona &: EngelbrecM 
(|1986I ) . They found periodic light curve variations 
in several Be stars in NGC 3766 (our Nos. 47, 
92, 133, 154, and 200). Additional Be stars, 
our Nos. 98, 130, 198 also exhibit light curve 
variations that are possibly periodic. In several 
cases, the shapes and possibly periods of the light 
curves change with timescales of only a few weeks. 
Balona et al.l (|l99l[ ) argue that this is due to mag- 
netic fields, but the light curves are not quali- 
tatively similar to variations observed in other 
magnetic stars, such as th e prototype a Ori E 
(jOksala fc Townsendl [2007h . The changing light 
curves among NGC 3766 members are probably 
more consistent with changes in active pulsation 
modes over short timescales. 

8. Conclusions 

Our spectroscopic analysis of NGC 3766 has re- 
vealed that Be stars may be much more common 
than we originally t hought. In our photorn etric 
study of NGC 3766 (iMcSwain fc Giejl2005bl) . we 
found up to 13 Be stars (5 definite, 8 uncertain) 
out of an expected 191 B-type stars, not count- 
ing the one Be star that saturated our photome- 
try. The new total of 16 Be stars is 23% greater. 
Among these 16 Be stars, 2-5 of them appear to 
have almost no disk at any given time, and an 
additional 2-4 have extremely subtle emission in 
their Ha line profile that could easily be mistaken 
for other phenomena (such as NRP manifesting 
themselves as bumps moving across the line or 
SB2 line blending). Therefore 25-50% of the Be 
stars may go undetected in a single spectroscopic 
observation, and photometric snapshots are even 
less likely to discern such weak emitters. We note 
four stars (Nos. 27, 45, 49, and 77) that were found 
to be possible or likely Be stars in the photometric 
study by Shobbrook (1985, 1987), but they never 
showed emission during our observations and thus 
remain unconfirmed. The existence of transitory, 
weak disks (especially Nos. 130 and 196) could 
mean that many more Be stars are waiting to be 
discovered. 

For our total sample of 48 Southern open clus- 



ters in our photom etric survey, we found a low Be 
fraction of 2 - 7% (jMcSwain & Giesll20n5bh . Con- 
sidering the very weak disks that are observed in 
NGC 3766 and the exceptionally high variability 
among the cluster's Be population, the total frac- 
tion of Be stars could be much greater. We are 
currently performing a similar spectroscopic study 
of several other clusters from our survey, and we 
will address those results in a future paper. 

While the Be stars of NGC 3766 are not dis- 
tinguishable from normal B-type stars by their 
evolutionary states, they do form a population 
of rapidly rotating stars. With two exceptions, 
their measured velocities are consistent with a uni- 
form population of rapid rotators having V = 
0.7 — 0.8 Vcrit- Gravitational darkening and weak 
emission in the He I lines may mean that these 
velocities are underest imated by as much as 33% 
(jTownsend et al.ll2004[ ). so the true T^ot is proba- 
bly at least 0.84 Vcrit- From the measured changes 
in the disks' masses and angular momenta, NRP 
are a capable source for the mass flow into the 
equatorial plane. The pulsations may be a tran- 
sitory phenomenon, however, and the variable na- 
ture of the Be stars probably reflects dramatic 
changes in the surface activity. 
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Fig. 1. — Ha profiles, labeled by MG ID num- 
ber. Spectra from 2003, where available, are 
shown with dotted lines, spectra from 2005 are 
shown with dot-dashed lines, spectra from 2006 are 
shown with dashed lines, and spectra from 2007 
are shown with solid lines. 
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Fig. 4. — Toff measured for the B-type tempera- 
ture calibratfon stars from this work (diamonds) 
and from Napiwotzki (199 3; triangles) compared 
to their calculated Teaiona ( Balonal[l984l ). A hnear 
fit to the two temperature scales (solid line) and 
the 1:1 agreement (dotted line) are also shown. 
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Fig. 6. — A contour plot of the la, 2(7, and 3a 
errors in the He I A4388 hne fit for No. 31. We 
adopted the center of the la error region as the 
best fit value, and the extent of this region indi- 
cates the sizes of the error bars for Tcff and \ogg. 



Fig. 5. — Stromgren cq index plotted against the 
(3 magnitude for the B-type temperature calibra- 
tion stars from this work (diamonds). Be stars 
are also plotted (filled diamonds) to illustrate the 
contamination in /3 due to their disk emission. We 
also show the ( cn, 0) relations for lurninosit y class 
V and III stars ( Balona fc Shobbrookl[l98i :' dotted 
lines) and our linear fit to this evolving population 
(solid line). 
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Fig. 7. — Toff and log.qp oiar are plotte d with 
the evolutionary tracks of Schaller et alj (Il992l) 
(solid lines) and isochrones of lLeieune fc Schaerer 
( 200l[ ) {dotted lines). The zero-age MS mass of 
each evolutionary track is labeled along the bot- 
tom, and the age of each isochrone is labeled 
along the top. Normal B-type stars are shown as 
open diamonds while Be stars are filled diamonds, 
and each symbol size is proportional to the star's 
V sini. 





Fig. 8. — Cumulative distribution function of 
V sini for the Be stars (dashed line) and the nor- 
mal B-type stars (dotted line). Theoretical dis- 
tributions of rapidly rotating stars with random 
orientation and V = 0.7 Vdit and 0.8 Vdit are also 
shown (solid lines; assumes the mean value for the 
Be stars, V^^it — 386 km s~^). 
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Table 1 
Journal of Spectroscopy 



UT 




Range 


Resolving Power 


Number of 




Telescope + 


Slit Plate 






Dates 




(A) 


(A/AA) 


Spectra 




Spectrograph 


{/xm) Grating 


Filter 


Detector 


2003 Mar 21 - 


- 22 


5490 - 6790 


1800 


20 


CTIO 


1.5ni + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2005 Feb 2 




4100 - 6900 


1900 


47 


OTIC 


Blanco 4m + Hydra 


200 KPGL3/1 




SITe 4Kx2K 


2006 May 13 




3790 - 4708 


3170 


38 


CTIO 


Blanco 4m + Hydra 


KPGLD/2 


BG39 


SITe 4Kx2K 


2006 May 14 


- 15 


5125 - 8000 


1560 


47 


CTIO 


Blanco 4iii -r Hydra 


KPGL3/1 




SITe 4Kx2K 


2007 May 4 - 


5 


5125 - 8000 


2000 


45 


CTIO 


Blanco 4m + Hydra 


200 KPGL3/1 




SITe 4Kx2K 


2007 Jan 20 




5650 - 6790 


1700 


10 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Jan 29 




5650 - 6790 


1700 


2 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Feb 2 




5650 - 6790 


1700 


8 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Apr 26 




5650 - 6790 


1700 


5 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Jun 9 




5650 - 6790 


1700 


1 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Jun 30 




5650 - 6790 


1700 


9 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Jul 3 




5650 - 6790 


1700 


2 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 


2007 Jul 27 - 


28 


5650 - 6790 


1700 


7 


CTIO 


1.5m + Cassegrain 


47/1 


GG495 


Loral IKxlK 



Table 2 
Ha Equivalent Widths 



MG 


WEBDA 


Wx (A) 


Wx (A) 


Wx (A) 


Wx (A) 


Wx (A) 


Wx (A) 


Wx (A) 


Wx (A) 


ID 


ID 


(2003 Mar) 


(2005 Feb) 


(2006 May) 


(2007 Jan/Feb) 


(2007 Apr) 


(2007 May) 


(2007 June) 


(2007 July) 


2 






3.24 








3.93 






1 a 

ID 


1 RQ 




A OS 


A Oft 






o.o4 






23 


■ ' ' 


' ' ' 


9.22 


9.37 






7.59 


' ' ' 




25 


291 


3.62 


3.18 


5.30 


— i.oo, — i.UD 


0.42 


1.04 


1.42 




27 


146 


5. 14 


5.38 


5.69 






4.74 






Oi 


1 

iOi 


O Q1 


3.99 


0.57 


Q nA OlA 
— 0.04, — o.y4 


R RA 


— 6.80 






36 


13 




8.46 


8.56 






6.93 






/1 1 


1 Qn 




4.88 


5.28 






/I /1 1 






42 


178 


' ' „ 


6.07 


6.02 






5.22 






45 


8 


4.69 


4.98 


5.10 






4.35 






/I? 


1 K 
ID 


1 A QQ 


1 1 QQ 


1 1 OQ 


ft 1 1 




R fil 
— D.oi 


— O.Oi 


— 4. DO 


/10 

4y 






A IT 


A QQ 






A OK 
4. ZD 






04 


125 




5.18 


5.61 






4.04 






55 






5.03 


5.12 






4.59 






57 


A 

24 




5.67 


5.86 






4.72 






61 


on 




4.09 


3.54 


O.OD 




3.56 


4.41 


1 on 


69 






2.89 


3.19 






2.81 




' " 


72 


4 


7.19 


8.27 


8.65 


7. 14 




7.08 




7.25 


73 


26 


3.65 


0.42 


3.28 


1.27 




— 0.35 




—0.92 


77 


195 


4.85 


4.61 


5.72 






4.73 






83 


27 


0.53 


4.07 


0.84 


2.92 


3.72 


2.93 






92 


1 


2.06 


2.31 


2.73 


3.32 




3.41 


2.66 




94 


194 




5.94 


5.91 






4.94 






96 


45 




6.39 


6.87 






6.26 






98 


36 


4.99 


3.10 


0.21 


4.43 


4.36 


4.14 




5.05 


101 


34 




7.40 


7.44 












118 






7.40 


7.03 






5.73 






119 


81 


4.76 


-5.48 


-3.19 


-6.53 




-8.67 




-12.83 


126 






7.46 


8.49 






6.86 






127 


53 


-7.21 


-4.95 


-4.00 


-4.19 




-5.47 


-5.60 


-12.81 


129 






8.66 


7.94 






6.84 






130 


67 




2.00 


5.04 


4.23 




4.44 


4.38 




133 


63 


4.40 


-1.15 


-9.14 


-11.62 




-12.72 


-12.40 




139 


204 


8.48 


5.46 


4.24 


4.51 


4.50 


4.42 


5.54 




154 


88 


-33.46 


-40.07 


-33.56 


-34.37 




-33.95 




-36.16 


155 






7.92 


7.90 






6.38 






161 


70 




3.92 


4.19 






3.61 






162 






9.15 


9.23 






7.45 






170 


94 




4.81 


5.26 






4.17 






173 


233 




7.79 


7.92 






6.62 






175 


218 




6.67 


6.74 






6.01 






178 


213 




6.32 


7.71 






6.35 






190 


260 




6.81 


6.86 












196 


239 


4.52 


4.32 


4.39 






1.46 






197 


253 




7.65 


8.85 


7.04 




7.09 




7.66 


198 


264 


-39.49 


-43.40 


-53.71 


-52.77 




-53.44 


-54.69 




200 


240 


-8.84 


-6.66 


-4.62 


-4.88 




-5.19 


-5.76 
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Table 3 

Physical Parameters of Cluster Members 



MG 


V sin i 


A.V sin i 


Toft 


ATeff 










AM* 


R* 






ID 


(km s"^) 


{km s-^) 


(K) 


(K) 


logs 


Alogg 


log gpolar 




(Me) 


(R@) 


(^0) 


(km s"i) 



2" 


307 


5 


17071 


450 


3.10 


0.05 


3.57 


7.0 


0.3 


7.2 


0.6 


352 


16" 


144 


7 


19420 


300 


3.48 


0.02 


3.62 


8.6 


0.2 


7.5 


0.5 


380 


27" 


287 


6 


16033 


200 


3.83 


0.02 


4.10 


4.9 


0.1 


3.3 


0.8 


436 


36" 


283 


22 


11931 


145 


4.02 


0.05 


4.27 


3.0 


0.0 


2.1 


0.1 


424 


41" 


84 


7 


17900 


300 


3.78 


0.03 


3.84 


6.6 


0.1 


5.1 


0.8 


406 


42=' 


239 


5 


16100 


100 


3.90 


0.02 


4.11 


4.9 


0.0 


3.2 


0.8 


438 


49" 


252 


9 


17400 


100 


3.75 


0.02 


3.99 


6.0 


0.1 


4.1 


0.9 


431 


54" 


192 


9 


17340 


250 


3.82 


0.02 


3.99 


5.9 


0.1 


4.1 


0.9 


430 


55" 


121 


6 


18000 


200 


3.78 


0.03 


3.88 


6.6 


0.1 


4.9 


0.8 


413 


57" 


238 


8 


16760 


250 


3.98 


0.03 


4.17 


5.1 


0.2 


3.1 


1.0 


459 


ei"-' 


331 


8 


18883 


350 


3.23 


0.02 


3.71 


7.7 


0.3 


6.4 


0.6 


390 


77" 


338 


10 


16000 


200 


3.80 


0.05 


4.15 


4.8 


0.1 


3.1 


0.9 


445 


94" 


177 


5 


15650 


150 


3.78 


0.03 


3.94 


5.0 


0.1 


4.0 


0.5 


401 


96" 


121 


5 


13937 


171 


4.31 


0.04 


4.37 


3.6 


0.1 


2.0 


0.3 


472 


101=* 


309 


14 


12753 


192 


3.84 


0.05 


4.14 


3.5 


0.1 


2.6 


0.4 


409 


118* 


242 


9 


13353 


191 


3.89 


0.05 


4.11 


3.7 


0.1 


2.8 


0.5 


409 


126* 


83 


9 


12399 


318 


3.79 


0.09 


3.87 


3.8 


0.1 


3.7 


0.4 


358 


129* 


97 


27 


12034 


203 


3.77 


0.06 


3.87 


3.7 


0.1 


3.7 


0.3 


355 


155" 


231 


11 


12694 


153 


4.05 


0.04 


4.22 


3.3 


0.1 


2.3 


0.5 


425 


161" 


73 


7 


18400 


250 


3.43 


0.02 


3.49 


8.3 


0.2 


8.5 


0.4 


351 


162" 


194 


11 


11513 


148 


3.99 


0.05 


4.12 


3.0 


0.1 


2.5 


0.2 


390 


170" 


65 


6 


18060 


250 


3.78 


0.03 


3.82 


6.8 


0.1 


5.3 


0.7 


403 


173" 


171 


7 


14210 


192 


4.21 


0.04 


4.30 


3.8 


0.1 


2.3 


0.1 


459 


175" 


252 


13 


13823 


193 


4.01 


0.05 


4.21 


3.8 


0.1 


2.5 


0.6 


435 


178" 


160 


12 


12504 


141 


3.92 


0.04 


4.05 


3.5 


0.1 


2.9 


0.3 


390 


197"''' 


198 


15 


12402 


182 


4.15 


0.06 


4.25 


3.2 


0.1 


2.2 


0.5 


427 



25° 


261 


6 


18995 


301 


4.02 


0.09 


4.02 


6.6 


0.1 


4.1 


0.9 


450 


31'* 


197 


5 


17834 


400 


3.82 


0.10 


3.99 


6.2 


0.2 


4.2 


1.0 


435 


47" 


190 


5 


18399 


301 


3.30 


0.09 


3.30 


9.3 


0.1 


11.3 


2.3 


323 


73'' 


296 


5 


18274 


301 


3.49 


0.09 


3.49 


8.2 


0.1 


8.5 


1.8 


349 


83" 






18817 


301 


3.31 


0.09 


3.31 


9.9 


0.1 


11.5 


2.4 


329 


92" 


214 


17 


18725 


301 


3.34 


0.09 


3.34 


9.5 


0.1 


10.8 


2.3 


332 


98" 






16890 


301 


3.84 


0.10 


3.84 


6.2 


0.1 


4.9 


1.0 


398 


119= 






17792 


301 


3.75 


0.09 


3.75 


6.8 


0.1 


5.8 


1.2 


387 


127'* 


347 


5 


17687 


550 


3.61 


0.12 


4.01 


6.1 


0.3 


4.0 


1.0 


437 


130° 


285 


7 


17519 


301 


3.69 


0.09 


3.69 


6.9 


0.1 


6.2 


1.3 


375 


133"^ 






18564 


301 


3.53 


0.09 


3.53 


8.3 


0.1 


8.1 


1.7 


358 


139° 


343 


12 


15945 


301 


3.95 


0.10 


3.95 


5.2 


0.1 


4.0 


0.8 


406 


154° 


112 


29 


13254 


1254 


3.29 


0.18 


3.29 


5.6 


0.3 


8.8 


1.8 


282 


196" 


165 


5 


19660 


250 


3.73 


0.02 


3.87 


7.6 


0.2 


5.3 


1.1 


426 


198'= 


251 


6 


19580 


301 


4.00 


0.09 


4.00 


7.0 


0.1 


4.4 


0.9 


449 


200'= 


236 


12 


16301 


301 


3.51 


0.10 


3.51 


6.8 


0.1 


7.5 


1.6 


337 



^Teff and log g measured from H7 line fit. 

^Parameters are unreliable since the star is a suspected SB2. 
*^Teff and log 5 measured from Stromgren photometry. 
*^Teff and log 5 measured from He I line fits. 
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Table 4 

Estimated Sizes of Be Star Disks 



MG 


HJD- 


Assumed i 


log Po 








ID 


2,450,000 


V/Vcrit (deg) 


(g cm-3) 


^disk/ R* 


(1O-"M0) 


(1O-"M0 yr-i) 


25 




0.7 64.9 


— 12.2 


2.3 


41.9 






3403.79 




-12.1 


2.5 


44.3 


1.3 




3870 




-12.3 


1.6 


30.6 


-10.7 




/lino c £2 A 




— 12.0 


3.8 


66. 1 


51.7 




/linn T /I Q 




— 12.0 


3.7 


63.7 


— 92.7 




4zib.4yo 




— 12.0 


3.3 


57.5 


—26.3 








— 12.1 


3.1 


54.7 


— 110.0 








— 12.1 


3.0 


53.0 


— 17.8 


31 


2720.744 


0.7 40.5 


-12.1 


2.5 


19.8 






3403.79 




— 12.2 


2.1 


17.2 


— 1.4 




3870 




— 12.0 


3.0 


24.6 


5.8 








— 11.9 


3.7 


31.4 


9.9 




/linn 1 




— 11.9 


4.1 


35.0 


142.4 




/I o 1 ii? rr 1 o 

4zib.oiz 




— 11.9 


4.1 


35.7 


3.3 




4225.66 




-11.9 


4.2 


35.9 


6.7 


47 


z7iy.o57 


0.7 40.8 


— 11.8 


4.5 


72.3 


' ■ ' 




o40o. * y 




— 11.8 


4.4 


70.3 


— 1.1 




3870 




— 11.8 


4.4 


70.2 


0.0 




/I 1 o o c ^ c 

4ioo.bob 




— 11.9 


4.2 


66.3 


—5.4 








— 11.9 


4.1 


64. 1 


—8.8 




/I no 1 /I o 




— 11.9 


4.0 


61.7 


— 15.7 




4309.466 




-11.9 


3.9 


59.5 


-28.5 


73 


z / iy.oo 1 


0.7 73.1 


— 12.2 


2.5 


16.7 


' ■ ■ 




o40o. * y 




— 12.0 


3.5 


23.0 


3.4 




3870 




— 12. 1 


2.6 


17.5 


—4.3 




,1 1 on ^ T T 

4iz0.7ii 




— 12.1 


3.3 


21.5 


5.8 




4225.66 




-12.0 


3.7 


24.3 


9.9 




4308.554 




-12.0 


3.9 


25.3 


4.1 


92 


OTT n ii? n n 

z719.62y 


n ^ Aon 

0.7 48.9 


— 12.1 


2.7 


50.4 


' ■ ■ 




o4Uo. /y 




— 12.1 


2.6 


49. 1 


—0.7 




3870 




— 12.1 


2.5 


47.0 


— 1.7 




4133.668 




-12.1 


2.3 


43.9 


-4.3 




4225.66 




-12.2 


2.3 


43.4 


-1.9 




4281.510 




-12.1 


2.5 


47.3 


25.9 


127 


2719.759 


0.8 80.1 


-11.9 


5.6 


30.7 






3403.79 




-11.9 


5.2 


28.5 


-1.2 




3870 




-11.9 


5.0 


27.4 


-0.9 




4133.680 




-11.9 


5.0 


27.6 


0.3 




4225.66 




-11.9 


5.3 


29.1 


5.8 




4281.524 




-11.9 


5.3 


29.2 


0.9 




4309.495 




-11.8 


6.1 


33.8 


59.9 


130 


3403.79 


0.7 63.5 


-12.1 


3.0 


11.9 






3870 




-12.3 


1.8 


7.7 


-3.3 




4120.750 




-12.2 


2.2 


9.0 


1.9 




4225.66 




-12.2 


2.1 


8.7 


-1.1 




4281.540 




-12.2 


2.1 


8.8 


0.6 


139 


2720.796 


0.9 50.6 




0.0 


0.0 






3403.79 




-12.3 


1.8 


1.7 


0.9 




3870 




-12.2 


2.3 


2.2 


0.4 




4120.770 




-12.2 


2.2 


2.1 


-0.1 
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Table 4 — Continued 



MG 


HJD- 


Assuincci ■/ 


log pv, 








ID 


2, loO.OUO 




(g cm ■'■) 




(10 "--U . ) 


(in yr ') 




4216.632 




-12.2 


2.2 


2.2 


0.0 




4225.66 




-12.2 


2.3 


2.2 


0.0 




4281.572 




-12.3 


1.8 


1.7 


-3.2 


196 


3870 


0.7 33.5 




0.0 


0.0 






4225.66 




-12.1 


2.6 


15.0 


15.4 


198 


2720.658 


0.7 69.9 


-11.9 


5.0 


253.4 






3403.79 




-11.9 


4.9 


248.9 


-2.4 




3870 




-11.9 


4.7 


238.7 


-8.0 




4133.732 




-11.9 


4.7 


239.5 


1.1 




4225.66 




-11.9 


4.7 


238.9 


-2.4 




4281.589 




-11.9 


4.7 


237.8 


-7.1 


200 


2720.637 


0.7 46.2 


-11.8 


4.8 


9.9 






3403.79 




-11.9 


4.6 


9.4 


-0.3 




3870 




-11.9 


4.3 


8.7 


-0.5 




4133.741 




-11.9 


4.3 


8.8 


0.1 




4225.66 




-11.9 


4.4 


8.9 


0.4 




4281.599 




-11.9 


4.5 


9.1 


1.2 
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